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Abstract. This paper considers the minimax regret 1-median problem 
in dynamic path networks. In our model, we are given a dynamic path 
network consisting of an undirected path with positive edge lengths, 
uniform positive edge capacity, and nonnegative vertex supplies. Here, 
each vertex supply is unknown but only an interval of supply is known. 
A particular assignment of supply to each vertex is called a scenario. 
Given a scenario s and a sink location a: in a dynamic path network, let 
us consider the evacuation time to r of a unit supply given on a vertex 
by s. The cost of x under s is defined as the sum of evacuation times to 
X for all supplies given by s, and the median under s is defined as a sink 
location which minimizes this cost. The regret for x under s is defined 
as the cost of x under s minus the cost of the median under s. Then, the 
problem is to find a sink location such that the maximum regret for all 
possible scenarios is minimized. We propose an 0(n®) time algorithm for 
the minimax regret 1-median problem in dynamic path networks with 
uniform capacity, where n is the number of vertices in the network. 

Keywords: minimax regret, sink location, dynamic flow, evacuation 
planning 


1 Introduction 

The Tohoku-Pacific Ocean Earthquake happened in Japan on March 11, 2011, 
and many people failed to evacuate and lost their lives due to severe attack by 
tsunamis. From the viewpoint of disaster prevention from city planning and evac¬ 
uation planning, it has now become extremely important to establish effective 
evacuation planning systems against large scale disasters in Japan. In particular, 
arrangements of tsunami evacuation buildings in large Japanese cities near the 
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coast has become an urgent issue. To determine appropriate tsunami evacuation 
buildings, we need to consider where evacuation buildings are located and how 
to partition a large area into small regions so that one evacuation building is 
designated in each region. This produces several theoretical issues to be consid¬ 
ered. Among them, this paper focuses on the location problem of the evacuation 
building assuming that we fix the region such that all evacuees in the region are 
planned to evacuate to this building. In this paper, we consider the simplest case 
for which the region consists of a single road. 

In order to represent the evacuation, we consider the dynamic setting in 
graph networks, which was first introduced by Ford et al. m- In a graph net¬ 
work under the dynamic setting, each vertex is given supply and each edge is 
given length and capacity which limits the rate of the flow into the edge per 
unit time. We call such networks under the dynamic setting dynamic networks. 
Unlike in static networks, the time required to move supply from one vertex to 
a sink can be increased due to congestion caused by the capacity constraints, 
which require supplies to wait at vertices until supplies preceding them have 
left. In this paper, we consider the flow on dynamic networks as continuous, 
that is, each input value is given as a real number, and supply, flow and time 
are defined continuously. Then each supply can be regarded as fluid, and edge 
capacity is defined as the maximum amount of supply which can enter an edge 
per unit time. The 1-sink location problem in dynamic networks is defined as the 
problem which requires to find the optimal location of a sink in a given dynamic 
network so that all supplies are sent to the sink as quickly as possible. 

In order to evaluate an evacuation, we can naturally consider two types of 
criteria: eompletion time criterion and total time criterion. In this paper we 
adopt the latter one (for the former one, refer to [12115117118) 1. We here define 
a unit as an infinitesimally small portion of supply. Given a sink location x in 
a dynamic network, let us consider an evacuation to x starting at time 0 and 
define the evacuation time of a unit to x as the time at which the unit reaches x 
in the evacuation. The total time for the evacuation to x is defined as the sum of 
evacuation times over all infinitesimal units to x. Then, the minimum total time 
for all possible evacuations to x could be the criterion for the optimality of sink 
location, which we adopt. Given a dynamic network, we define the 1-median 
problem as the problem which requires to find a sink location minimizing the 
minimum total time, and the optimal solution is called the median. 

Although the above criterion is reasonable for the sink location, it may not 
be practical since the number of evacuees in an area may vary depending on 
the time (e.g., in an office area in a big city, there are many people during the 
daytime on weekdays while there are much less people on weekends or during 
the night time). So, in order to take into account the uncertainty of population 
distribution, we consider the maximum regret for a sink location as another eval¬ 
uation criterion assuming that for each vertex, we only know an interval of vertex 
supply. Then, the minimax regret 1-median problem in dynamic path networks is 
formulated as follows. A particular assignment of supply to each vertex is called 
a scenario. Here, for a sink location x and a scenario s, we denote the minimum 
total time by <P^{x). Also let m® denote the median under s. The problem can 
be understood as a 2-person Stackelberg game as follows. The first player picks 




a sink location x and the second player chooses a scenario s that maximizes the 
regret defined as <P^{x) — The objective of the first player is to choose x 

that minimizes the maximum regret. 

Related to the minimax regret facility location in graph networks, especially 
for trees, some efficient algorithms have been presented by |2|3|5|6 |7|9j. For dy¬ 
namic networks, Cheng et al. [5] first studied the minimax regret 1-eenter problem 
in path networks, which requires to find a sink location in a path that minimizes 
the maximum regret where the completion time criterion is adopted instead of 
the total time one. They presented an O(nlog^n) time algorithm. Higashikawa 
et al. m improved the time bound by [5] to 0(n log n), and also Wang [TO] 
independently achieved the same time bound of O(nlogn) with better space 
complexity. Very recently, Bhattacharya et al. [1] have improved the time bound 
to 0{n). The above problem was extended to the multiple sink location version 
by Arumugam et al. [1]. For the minimax regret fc-center problem in dynamic 
path networks with uniform capacity, they presented an 0{kn^ logn) time algo¬ 
rithm, and this time bound was improved to 0{kn^) recently |12j . On the other 
hand, for dynamic tree networks, only the minimax regret 1 -center problem was 
solved in 0{n^ log^ n) time |14I16) . 

This paper first considers the minimax regret median problem in dynamic 
networks while all the above works for dynamic networks treated center prob¬ 
lems. In this paper, we address the minimax regret 1-median problem in dynamic 
path networks with uniform capacity and present an 0(n^) time algorithm. 


2 Preliminaries 

2.1 Dynamic path networks under uncertain supplies 

Let P = {V, E) be an undirected path with ordered vertices V = {vi,V 2 , ■. ■, Vn\ 
and edges E = {ei, 62 ,..., e„_i} where = (vi,Vi+i) for i £ {1,..., n — 1}. Let 
Af = (P, 1,w,c,t) he a. dynamic network with the underlying path graph P; Hs a 
function that associates each edge with positive length k, w is a function that 
associates each vertex Vi with positive weight Wi, amount of supply at c is the 
capacity, a positive constant representing the amount of supply which can enter 
an edge per unit time; r is a positive constant representing the time required 
for a flow to travel a unit distance. In our model, instead of the weight function 
w on vertices, we are given the weight interval function W that associates each 
vertex Vi £ V with an interval of supply Wi = [w ~, wf ] with 0 < < wf. 

We call such a network M = (P, / , IF, c, r) with path structures a dynamic path 
network under uncertain supplies. 

In the following, we write p £ P to indicate that a point is a vertex of P or 
lies on one of the edges of P. For any point p £ P,we abuse this notation by also 
letting p denote the distance from vi to p. Informally we can regard P as being 
embedded on a real line with vi = 0. For two points p,q £ P satisfying p < q, 
let [p,cj\ (resp. [p,q), ip,q] and (p,q)) denote an interval in P consisting of all 
points X £ P such that p < x < q (resp. p<x<q, p<x<q and p < x < q). 






2.2 Scenarios 


Let S denote the Cartesian product of all Wi for i e {1,..., n}: 

n 

s = l[w.. ( 1 ) 

i=l 

An element of S, i.e., a particular assignment of weight to each vertex, is called 
a scenario. Given a scenario s G S, we denote by wf the weight of a vertex Vi 
under s. 


2.3 Total evacuation time 


In our model, the supply is defined continuously. We define a unit as an infinites¬ 
imally small portion of supply. Given a sink location x G P and a scenario s G 5, 
without loss of generality, an evacuation to x under s is assumed to satisfy the 
following assumptions. When a unit arrives at a vertex v on its way to x, it has 
to wait for the departure if there are already some units waiting for leaving v. 
All units waiting at v for leaving v are processed in the first-come first-served 
manner. 

For a given x G P and s G 5, let us consider an evacuation to x under s 
starting at time 0 and define the evacuation time of a unit to x under s as the 
time at which the unit reaches x. Let <?®(x) denote the sum of evacuation times 
over all infinitesimal units to x under s. Also let ^‘l{x) (resp. ^|j(x)) denote the 
sum of evacuation times to x under s for all units on [ui, x) (resp. (x, «„]). Then, 
<?®(x) is obviously the sum of ^^^(x) and ^|{(x), i.e., 

^^ix)=^Ux)+^U^). ( 2 ) 

Without loss of generality, we assume = 0 and = 0. 

We will show the formula of (x) that has been proved in |15I17] . Suppose 
that X is located in an open interval {vh,Vh+i) with 1 < < n — 1, i.e., x G eh- 

We here show only the formula of ^%{x) (the case of ^|i(x) is symmetric). First, 
let us define the vertex indices pi,..., pk inductively as 


Pi = argmax 



i) + 






3 G {pi-i -I- 1,..., /i} 


(3) 


where po = 0. Obviously pk = h holds. We then call a set of all units on 
[up;_j+i,UpJ i-th cluster, and a vertex Vp^ the head of i-th cluster (see Fig. [1]). 
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Fig. 1. Illustration of left clusters for x where i-th cluster is headed by a vertex 


















Also, for i S fc}, we define Oi as Gi = ^ which is called 

the weight of j-th cluster. The interpretation can be derived from mm as 
follows. The first unit on each Vp. does not encounter any congestion on its way 
to X. Here, although pi may not be uniquely determined by ([3]), we choose the 
maximum index as pi . By this assumption and ©, the following inequality holds 
for j G {2 ,..., k}\ 


T{Vp,-Vp,_^)>'^. (4) 

In other words, even if we transform the input so that all units on [vp^_^+i,Vp^) 
are moved to for i G {!,..., fc}, the sum of evacuation times over [vi,x) 
never changes. In the following, we call such a transformation the left-clustering 
for X and clusters obtained by the left-clustering for x left clusters for x (the 
right-clustering and right clusters are symmetric). Then as in |15I17) . 'P‘l{x) is 
written as 




(5) 


2.4 Minim 2 Lx regret formulation 

For a scenario s G <S, let be a point in P that minimizes ^‘^{x) over x G P, 
called the median under s. We now define the regret for x under s as 

=r(a:)-r(m"). (6) 

Moreover, we also define the maximum regret for x as 

^max(^ ) = max{i?®(x) I s G 5}. (7) 

If s = argmax{i?®(a:) | s G 5}, we call s the worst case scenario for x. The 
goal is to find a point x* G P, called the minimax regret median, that minimizes 
Rmax{x) over X G P, i.e., the objective is to 

minimize {i?max(a^) | x G P}. (8) 

2.5 Known properties for the fixed scenario case 

We here show some properties on the 1-median problem in a dynamic path 
network M = {P = (V, E), I, ir®, c, r) when a scenario s G 5 is given, which were 
basically presented in |15I17) . We first introduce the following two lemmas. 

Lemma 1 |15I17] For a scenario s G S, m® is at a vertex in V. 

Lemma 2 |15I17] For a scenario s G S, all ^®(ui) over i G {!,...,u} can be 
computed in 0{n) time in total. 


We then can see a corollary of these lemmas. 












Corollary 1 |15I17) For a scenario s G S, and can be computed in 

0(n) time. 

Now let us look at the formula of ©■ Even if x is moving on an edge Ci 
(not including endpoints Vi and Ui+i), the formation of left clusters for x does 
not change over x G Ci. Therefore, ^£(x) is a linear function of a; G Ci, and 
symmetrically, is also a linear function. For i G {1,... ,n — 1}, letting a® 

and bf be the values such that for x G et, 

<?®(a;) = a®x + 6®, (9) 

we can derive the following lemma from [T5[T7] . 

Lemma 3 For a scenario s G S, all a® and &® over i G {1, ... ,n — 1} can be 
computed in 0{n) time in total. 


3 Properties of worst case scenarios 

In this section, we show the important properties which worst case scenarios 
have. In our problem, a main difficulty lies in evaluating R^(x) over s G S to 
compute Rmaxix) even for a fixed x since the size of S is infinite. We thus aim 
to find a scenario set with a finite size (in particular, a polynomial size) which 
includes a worst case scenario for any x G P. In order to do this, we introduce 
a new concept, the gap between two points x,y G P under a scenario s G S, 
defined by 


T®(x,j/)=<l>®(x)-<Z>®(y). 

By Lemma [Hand the definition of ([5]), we have 

i?®(a;) = max{T®(a;,y) \ y gV}, 


( 10 ) 


( 11 ) 


and by d?]) and (fTTl) . 

l?max(a;) = max{max{T®(a;,?/) | y S E} | s e 5} 

= max{max{r'®(a;, y) | s £ <S} | y S V}. (12) 

From if we can compute max{F®(a;,y) | s £ 5} for a fixed pair {x,y) G 
P X y, Rmax(x) can also be computed by repeating the same maximization over 

y G V. We call a scenario that maximizes F^(x, y) for a fixed {x, y) a worst case 

scenario for {x,y). In the following, we show a scenario set of size 0{n) that 
includes a worst case scenario for a fixed {x,y), which implies a scenario set of 
size 0{n^) that includes a worst case scenario for a fixed x. 

3.1 Bipartite scenario 

We first introduce the concept of the bipartite scenario, which was originally 
introduced as the dominant scenario in msi- Let us consider a scenario s G S. 








A scenario s is said to be left-bipartite (resp. right-bipartite) if w® = w'j' (resp. 
w~) over j G {1,..., i} and re® = w~ (resp. ) over j S {z + 1,..., n} for some 
i G 1}. Obviously the number of such scenarios is 0{n). The authors 

of |8ll3j treated the minimax regret 1-center problem in dynamic path networks, 
which requires to find a sink location in a path that minimizes the maximum 
regret similarly defined as © where the completion time criterion is adopted 
instead of the total time one. They proved that for any point in an input path, 
at least one worst case scenario is left-bipartite or right-bipartite. 


3.2 Pseudo-bipartite scenario 

We here introduce the concept of the pseudo-bipartite scenario. A scenario s is 
said to be left-pseudo-bipartite (resp. right-pseudo-bipartite) if w® = w'j' (resp. 
w~) over j G {l,...,z — 1} and rc® = w~ (resp. ) over j G {z -I- 1,...,n} 
for some z G {2,...,n — 1}. In this definition, we do not care about the weight 
of a vertex Vi, called the intermediate vertex. Given a pseudo-bipartite scenario 
with the intermediate vertex we call intervals [z;i, Vi) and {vi, Vn] the left part 
and the right part, respectively. Let Sl (resp. Sr) denote a set of all left-pseudo- 
bipartite scenarios (resp. right-pseudo-bipartite scenarios). We then prove the 
following lemma (the proof is given in Appendix A). 

Lemma 4 Given a pair {x,y) G P xV satisfying y < x (resp. x < y), there 
exists a worst case scenario for (x, y) belonging to Sr (resp. Sr) such that y 
(resp. x) is in the left part and x (resp. y) is in the right part. 


3.3 Critical pseudo-bipartite scenario 

By LemmalU we studied the property of a worst case scenario for a fixed (x, y) G 
P X V, however the sizes of Sl and Sr are still infinite since the weight of the 
intermediate vertex in a pseudo-bipartite scenario is not fixed. In the rest of this 
section, we focus on the weight of the intermediate vertex in a pseudo-bipartite 
scenario which is worst for (x,y). 

Given a pair (x, y) G PxV satisfying y < x, let us consider a scenario s G Sr 
such that the intermediate vertex is Vi and y < Vi < x. Suppose that the weight 
of Vi is set as the minimum, i.e., wf = w~. Performing the right-clustering for y 
under s (mentioned in Section 12.31) , we will get right clusters for y such that for 
I G {l,...,/c}, the head of Tth cluster is pi and the weight of l-th cluster is ai. 
Then, suppose that the intermediate vertex Vi belongs to j-th cluster. 

Now let us increase the weight of Vi, little by little, without changing the 
weight of any other vertex. Let s(w) be a scenario in Sr such that the interme¬ 
diate vertex is Vi whose weight is w G [w~ , wf ]. Suppose that when the weight 
of Vi reaches some value ui, the following equality holds: 


TiVp._, -Vpft 


<7j P{UJ-W^ ) 

c 


(13) 


Note that (Jj — ) corresponds to the weight of j-th cluster under s(a;). At 

that moment, referring to (|1]), (j —l)-th cluster is merged to j-th cluster. We then 





call s{uj) a critical left-pseudo-bipartite scenario for y. Also, s{w~) and s{wi^) 
are assumed to be critical left-pseudo-bipartite scenarios for y even if any merge 
does not occur at those moments. Critical right-pseudo-bipartite scenarios for y 
are symmetrically defined. Let Sy denote a set of all critical left-pseudo-bipartite 
scenarios and critical right-pseudo-bipartite scenarios for y, and S* = [Jy^ySy. 
We will show two lemmas (the proof of Lemma [5] is given in Appendix B). 

Lemma 5 Given a pair {x, y) € P x V, there exists a worst case scenario for 
{x,y) belonging to Sy. 

Lemma 6 For a vertex y €V, the size of Sy is 0{n), and all scenarios in Sy 
can be computed in 0(n) time. 

Proof. We first prove that the number of critical left-pseudo-bipartite scenarios 
for y is 0{n) (the case of critical right-pseudo-bipartite scenarios is symmetric). 
Suppose that y = vj. For i G {j -|- 1,...,n} and w G [rc”, w)*"], let s{i,w) be a 
scenario in Sl such that the intermediate vertex is Vi whose weight is w. Here, 
let us define the order between two scenarios s{i,w) and s(i',w'): s{i,w) -< 
s(j', w') holds if and only if (a) i < i' or (b) i = i' and w < w'. For i G 
{j -\- l,...,n}, we also define pi and qi as follows. Let pt be the number of 
critical left-pseudo-bipartite scenarios for y such that the intermediate vertex is 
Vi (including s{i,w~) and s{i,wf)). Let qi be, under a scenario s{i,wf), the 
number of right clusters for y that follow a cluster including Vi. 

Let us consider computing all critical left-pseudo-bipartite scenarios for y in 
ascending order, and suppose that the weight of Vi now increases from w~ to 
wil. While it increases, since pi — 2 critical left-pseudo-bipartite scenarios for 
y occur (except s{i,w~) and s{i,wf)), and at each such scenario, one or more 
clusters are merged into one of Vi, at least Pi — 2 clusters are merged into one of 
Vi in total. We thus have qi < qi-i — {pi — 2) for i G {j -\-1,..., n}, i.e.. 

Pi < qi-i - (?i -b 2. (14) 

Note that the total number of critical left-pseudo-bipartite scenarios for y is 
exactly 1 -b Y17=j+APi “ !)• (O, we have 


X] (P* “ 1) ^ X! + 1) = Qj -qn + in- j), (15) 

i=j + l 

which is 0{n) since qj < n — j and qn = 0. 

In the rest of the proof, we show that all critical left-pseudo-bipartite scenar¬ 
ios for y = Vj can be computed in 0{n) time. Recall that all critical left-pseudo- 
bipartite scenarios for y are computed in ascending order. The algorithm first 
gets s(j -b 1, and performs the right clustering for y under s{j -b 1, 

As claimed in msm, it is easy to see that the right clustering for a fixed y can 
be done in 0(n) time. 

Suppose that for particular i G {j -\- 1,.. . ,n} and to G [wf ,wf], s{i,uj) is 
critical for y, and the algorithm has already obtained s{i, oj) and the right clusters 
for y. We then show how to compute the subsequent critical left-pseudo-bipartite 




scenario. Let Cy be a right cluster for y including Vi and Cy be a right cluster for 
y immediately following Cy. Also, let py (resp. p'y) be the index of a vertex that 
corresponds to the head of Cy (resp. c'y), and Oy (resp. a'y) be the weight of Cy 
(resp. c'y). 

There are two cases: [Case 1] w < wf ; [Case 2] uj = wf. For Case 2, we notice 
that s{i + is equivalent to s{i,w^). Therefore, this case immediately 

results in Case 1 by letting i be i + 1 and oj be (although a right cluster 
for y including Vi+i may be c^, not Cy). We thus consider only Case 1 in the 
following. 

The algorithm will compute the subsequent critical left-pseudo-bipartite sce¬ 
nario s{i,uj') where uj' satisfies to < uj' < wf. In order to compute w', the 
algorithm test if there exists w G such that 


"Kl -^Py) = 


ay + [w- w) 


(16) 


which is similar to (nsi). If yes, for such w, the algorithm returns ui' = w and 
updates the right clusters for y by merging c'y into Cy. Otherwise, uj' = wf 
is just returned. Such testing and updating are done in 0(1) time. Since the 
number of critical left-pseudo-bipartite scenarios for y is 0(n) and each of those 
is computed in 0(1) time, we completes the proof. □ 


By ((HI), we have a corollary of Lemma |S] 

Corollary 2 Given a point x € P, there exists a worst case scenario for x 
belonging to S*. 

Also, a corollary of Lemma [5] immediately follows. 

Corollary 3 The size of S* is O(n^), and all scenarios in S* can be computed 
in O(n^) time. 


4 Algorithm 

In this section, we show an algorithm that computes the minimax regret median, 
which minimizes Rmaxix) over x G P. The algorithm basically consists of two 
phases: 

[Phase 1] Compute i?max(fz) over * G {1,..., n}, and 

[Phase 2] Compute min{i?max(a^) j a: G e^} over j G {1,..., n — 1}. 

After these, the algorithm evaluates all the 2n — 1 values obtained and finds the 
minimax regret median in 0(n) time. 

By Corollary[2j we only have to consider scenarios in S* to compute Rmaxix) 
for any x G P. Therefore, the algorithm computes all scenarios in S* in advance, 
which can be done in Oln^) time by Corollary [31 Subsequently, it computes all 
the values <?®(m'*) over s G S* for Phase 1 and Phase 2. By Corollaries [T] and [31 
this can be done in 0{n^) time in total. 

First let us see details in Phase 1. For a fixed scenario s G 5*, since all 
<T’^{vi) over i G {l,...,n} can be computed in 0(n) time by Lemma [3] and 



has already been computed before Phase 1, all R^{vi) over i G {1,..., n} 
can also be computed in 0{n) time (refer to (jH]))- After the algorithm obtains 
(vi),..., (vn) over s G 5* in 0(n^) time, for each i G n}, R^{vi) 

over s £ S* are evaluated to obtain Rmax{vi)- Thus, it is easy to see that Phase 

1 can be done in 0{n^) time in total. 

We next focus on Phase 2. As mentioned at the end of Section 1^751 for a fixed 
scenario s G 5*, <P^{x) is a linear function of x G for each i G {l,...,n — 1} 
(not including Vi and Xi+i). Therefore, i?®(x) is also linear for x G for each i. 

Referring to dH), a function i?'*(x) on an edge is written as 

i?"(x) =a|x + &|-<?*(m"). (17) 

Recall that ^^(m®) has already been computed. Then, by Lemma[21 i?®(x) on a 
over i G {1,..., n — 1} can be computed in 0{n) time. After the algorithm does 
the same computation over s G 5* in 0{n^) time, on each edge e^, we have 0{n?) 
linear functions i?®(x) over s £ S*. By the definition of o, min{i?i„ax(a;) | x G 
Ci} can be obtained by solving a linear programming problem in two dimensions 
with Oiji?) constraints, i.e., 

minimize y 

subject to a^x + bl — < y, Vs G <S* 

Vi < X < Vi+i. 

This problem can be solved in O(n^) time by [TU]. Repeating the same operations 
over * G {1,..., n — 1}, Phase 2 is completed in 0{n^) time. 

Theorem 1 The minimax regret 1-median problem in dynamic path networks 
with uniform capacity can be solved in 0{n^) time. 

5 Conclusion 

In this paper, we address the minimax regret 1-median problem in dynamic 
path networks with uniform capacity and present an O(n^) time algorithm. 
Additionally, this is the first work that treats the minimax regret facility location 
problem in dynamic networks where the total time criterion is adopted. Two 
natural questions immediately follow. The first one is whether we can reduce 
the number of scenarios to be considered. The other one is whether we can 
extend the problem to the fc-median version with fc > 2, or the problem in more 
general networks. 
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Appendices 

Appendix A: Proof of Lemma 


We only treat the case of {x,y) € P x V satisfying y < x since the other case is 
symmetric. 

We first prove the following claim. 

Claim 1 Given a pair {x,y) € P xV satisfying y < x, there exists a worst case 
scenario for {x,y) such that the weight of every vertex Vi G [ni,y] is w^. 

Proof of Claim{^ We here assume y > vi: ii y = vi, the proof is straightforward. 
Let Si be a worst case scenario for {x,y). If there are more than one worst case 
scenario, we choose the one such that all weights are lexicographically maximized 
in the order of ascending indices among all worst case scenarios. We prove by 
contradiction: suppose < w'l for some vertex Vi G [vi^y\. Under a scenario 
si, let us perform the left-clustering for x and y, respectively (refer to Section 
[^ . Performing the left-clustering for x, let Cx be a left cluster for x including 
Vi, Px be the index of a vertex that corresponds to the head of Cx, and Ux be 
the weight of Cx- Also, performing the left-clustering for y, Cy, py and Uy are 
similarly defined. When Cy evacuates to cc, by the definition of a cluster, the first 
unit of Cy does not encounter any congestion in the interval (vp^,y), however, 
it may encounter congestions at some vertices in [y-,x). If it does not encounter 
any congestion in (vp^,x), Cy is never merged to any other cluster for x, which 
implies px = Py and ax = ay. Otherwise, Cy is eventually merged into a left 
cluster for x headed by a vertex Vp^ € [y, x), i.e., Cx- We thus consider two cases: 
[Case I] Px = Py and ax = ery', [Case 2] py < y < px and ax > (Xy. 

Now let S 2 be a scenario obtained from si by increasing the weight of vt by 
infinitesimally small <5 > 0, i.e., = icU -|- S and Wj^ = Wj^ for j ^ i. We then 

show the following claim for the proof of Claim [T] 

Claim 2 While si changes to S 2 , the formation of left clusters for x (resp. y) 
remains the same. 


Proof of Claim\^ Performing the left-clustering for x under si, let be a left 
cluster for x immediately following Cx, and Px be the index of a vertex that 
corresponds to the head of c(,. Referring to dH), the following inequality holds: 


( ^ ^ 

"K. — = 


Si 


Then, for a sufficiently small <5 > 0, we have 


(18) 


- ''^Pi) > 


(E 


px 


^r)+<5 
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E Px 
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(19) 


The inequality of (1181) means that after performing the left-clustering for x under 
si, the first unit of c(, does not catch up with the last unit of Cx at Vp ^, and by 
HU), this remark also holds even for S 2 . Under S 2 , if c(, is not merged to Cx, 





any other merge never occurs. Thus, the formation of left clusters for x does not 
change, and similarly, it does not change for y. □ 


By Claim [2] and the definitions of ([2]) and ([5]) , we have 
(x) - (x) = (x) - (cc) 



( 20 ) 


and similarly. 


= ST{y-Vp^) + ^ + 


( 21 ) 


Also by the definition of OT, we have 


(x, y) - {x, y) = (x) - (y) - (<P*^ (x) - (y)) 

= (x) - (x) - (<?*= (y) - (y)). (22) 


From ((^ . (1^ and (1^ . we can derive 



(23) 


If Case 1 occurs, we can immediately see that the right side of (USD is greater 
than zero, i.e., C®2(x,y) > r^^{x,y), which contradicts that si is a worst case 
scenario for (x,y). 

If Case 2 occurs, performing the left-clustering for x merges Cy into Cx as men¬ 
tioned above, and then, all units on (vp ^, Vp^] are also merged into Cx- Therefore, 
if we consider an input such that all units on [vp ^, Vp ^) are moved to Vp ^, the 
first unit of Cy must catch up with the last unit of supply at i.e.. 


Since Cx includes Cy and all units on {vp^^VpJ^, we have 



(25) 


I - py-\-l 


Note that in (EH), the left side is less than the right side when Cx also includes 
some clusters for y following Cy. From (1^ . (|^ and (1^ . we can derive 


(x, y) - (x, y) > 6t{x - y) > 0, 


(26) 





which contradicts that si is a worst case scenario for {x,y). □ 

If we consider a worst case scenario for (a:, y) such that the weight of every 
vertex Vi G [ui,y] is and weights of all other vertices in (y, u„] are lexico¬ 
graphically minimized in the order of descending indices, the following claim is 
also proved in a similar manner as in the proof of Claim [TJ 

Claim 3 Given a pair {x,y) G P xV satisfying y < x, there exists a worst case 
scenario for {x, y) such that the weight of every vertex Vi G [ui, y] is and the 
weight of every vertex Vi G [x,Vn\ is w~. 

Now, let S3 be a worst case scenario for {x, y) such that the weight of every 
vertex Vi G [vi,y\ is the weight of every vertex Vi G [x, Vn] is w~, and weights 
of all other vertices in the open interval (y, x) are lexicographically maximized 
in the order of ascending indices. Then, S 3 G Sl can be proved. We prove by 
contradiction: there exist two vertices Vi,Vj G (y,x) satisfying i < j such that 
< wil and > w~. Let S4 be a scenario obtained from S3 by increasing 
the weight of Vi by infinitesimally small <5 > 0 and decreasing the weight of Vj 
by the same 5, i.e., = wP -\- 5, — S and wp = wp for k i,j. 

Then, we immediately see <P‘‘^{x) > <P^^{x) and ^^^(y) < <?®^(y), therefore 

(x) — (y) > (x) — (y), i.e., 

r^px,y)>r^px,y). ( 27 ) 

The inequality of (EH) implies that S 4 is also a worst case scenario for (x, y), which 
contradicts the lexicographical maximality of weights on the open interval (y, x) 
under S3. □ 


Appendix B: Proof of Lemma [5] 

For a fixed pair {x,y) G P x V satisfying y < x, let us consider a worst case 
scenario in iSl such that the intermediate vertex is Vi {y < Vi < x). We now 
consider the weight of Vi as a variable w G [w~, wi^], and let s{w) be a scenario 
in Sl such that the intermediate vertex is Vi whose weight is w. Then, let r(w) 
denote the gap between x and y under s{w), i.e., 

r{w) = - ^>"(“’Hy). ( 28 ) 

Suppose that s(w) is critically left-pseudo-bipartite for y when w = tvi, ... ,ujp, 
where p is a positive integer and w~ = uJi < ... < cOp = wf . In the following, 
we prove that a function r{w) is convex and piecewise-linear for w G [ujj,ujj+i] 
for every j G {1,... ,p — I}. 

Under a scenario s{w), let us perform the left-clustering for x and the right¬ 
clustering for y, respectively. Performing the left-clustering for x, let Cx{w) be a 
left cluster for x including Px{w) be the index of a vertex that corresponds 
to the head of Cx{w), and ax(w) be the weight of Cx. Also, performing the right¬ 
clustering for y, Cy{w), Py{w) and (Jy(w) are similarly defined. 

We first show the following claim. 


Claim 4 r{w) is continuous for w € 


Proof of Claim The statement is equivalent to 


lim r(w + (5) = riw) 

5^-+o 


Vw e [Wj 


lim rlw 

5^+0 


i5) = r{w) Vw G {w^ 


and 


(29) 

(30) 


We here prove (1^ (the case of (1501) is similarly treated). By (1^ . we only have 
to prove 


lim = r(“)(x) \/w G [w-,wt). (31) 

Note that, similarly as in Claim[5J while s{w) changes to s(i(; + (5), the formation 
of left clusters for x remains the same. Therefore, by the definitions of ([2]) and 
®, we have 

+ St{x - (32) 

which leads m by letting <5 go to positive zero. □ 


For an integer j G {l,...,p—1}, we consider the right-derivative of Piw) for 
w G [ujj,ujj+i), i.e., 


r+ {w) = lim 


r{w -I- (5) — T(w) 


(33) 


<5^+0 S 

Similarly to we have 

(y) =<!>«(“) (y) + 6t{v^^ („) -y) + ^. (34) 

From ((Mil . (IMll and (l34ll . we derive 


r{w + S)- r{w) = - ^>*(“+'5)(y) - |(|)«(“')(a;) - 


t- / , S ^rrxiw) — ay(w)} 

= 6t{x + y- Vp^^y,) - Vp^(y,)) -f- - -(35) 


and by ([551) and (155]) . 

^/ / N z X Gxiw) — (7y(w) 

r+{w) =T{x + y- Vp^i^y,) - Vp^i^y,)) P --. (36) 

We here notice that as w increases, Px(w) and Pyiw) never change (even if left 
clusters for x following Cx{w) and right clusters for y following Cy(w) are merged 
to Cx{w) and Cy{w), respectively). Also, since both of Cx{w) and Cy(w) include 







Vi, Oxiw) and (Jy{w) can be represented as follows: 


(37) 

(38) 


(Jx{w) = w + cjI.{w), and 
ay{w) = w + (J^y{w), 

where uKw) and cr^{w) are functions of w. In addition, aKw) increases only if a 
left cluster for x following Cx{w) is merged to Cx(w), i.e., cr^(rc) is an increasing 
step function, and cf^iw) is a constant function of w G [ujj,ujj+i) since the 
formation of right clusters for y does not change over w G [wj , w^+i) (recall the 
definition of critical left-pseudo-bipartite scenarios for y in Section r3.3|) . From the 
above observations, and (Ei), (EZl) and (I55)) . we derive that for w G w^+i), 

r'{w) = + const., (39) 

c 

which is an increasing step function. By this fact and the continuity of r{w) by 
Claim m we have the following claim. 

Claim 5 For an integer j G {1, ... ,p— 1}, r{w) is convex and piecewise-linear 
for w G [ujj,ujj+i]. 

By Claim 0 a solution that maximizes r(w) must be in {wi,..., ujp}, i.e., a 
worst case scenario for (x, y) is critically left-pseudo-bipartite for y. □ 



